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1.  INTRODUCTION 


Wall  pressure  fluctuations  produced  by  turbulent  boundary  layer  flows 
can  produce  structural  vibrations  and  noise.  In  order  to  predict  surface 
motion  details,  as  part  of  the  design  process,  realistic  wall  pressure  forc¬ 
ing  function  data  are  required  which  duplicate  the  statistical  properties  of 
turbulent  flows.  Ultimately,  a  coupled  fluid-structural  analysis  package 
which  employs  accurate  structural  modelling  and  the  Navier-Stokes  equations 
in  three-dimensions  is  desired  (Refs.  1,  2),  but  that  approach  presently  has 
limited  resolution  and  is  very  expensive.  Experimental  measurements  of 
P(x,z,t)  over  a  rectangular  region,  are  the  next  best  source  of  forcing 
function  data  for  the  design  process,  but  experimental  data  suffer  from 
limited  resolution,  expense  and  the  need  for  arbitrarily  resolvable  measure¬ 
ments  over  a  wide  range  of  geometrical  and  flow  conditions.  A  less  expen¬ 
sive  alternative  is  numerical  simulation.  A  Monte  Carlo  wall  pressure  simu¬ 
lation  has  been  developed  and  will  be  described  here. 

While  nonlinear  structural  systems  may  require  real  time  forcing  func¬ 
tion  input,  wall  pressure  data  transformed  into  frequency-wave  number  space 
is  preferred  for  most  linear  systems,  (Refs.  3-7).  Even  then  P(x,z,t)  data 
are  required  to  develop  appropriate  transformation.  Those  data  are  not 
available  in  general  and  most  structural  dynami cists  and  acousticians  have 
resorted  to  some  type  of  statistical  decomposition  assumption  which  permits 
prediction  of  surface  amplitude  spectra  but  may  not  yield  reliable  wave  num¬ 
ber  information.  Attempts  to  incorporate  more  realistic  turbulence  models 
have  been  reported  (Refs.  3,  7,  8),  but  uncertainty  still  exists  concerning 
interpretation  of  wall  pressure  statistics  (Refs.  6,  9)  and  there  may  be  a 
limit  to  the  validity  of  statistical  data  using  experimental  instruments. 

Questions  concerned  with  the  low  frequency  end  of  the  pressure  fluctua- 


tion  spectrum  and  limiting  local  root-mean  square  pressure  levels  remain 
(Ref.  9),  but  there  is  little  doubt  about  the  generic  shape  of  the  pressure 
spectrum  in  the  primary  energy  region  or  about  the  decreasing  correlation 
peak  with  separation  distance  between  spatial  locations.  As  will  be  shown 
subsequently,  the  statistics  produced  from  the  synthetic  pressure  distribu¬ 
tion  are  remarkably  similar  to  turbulent  wall  pressure  data.  While  the 
synthetic  case  does  not  replicate  the  physics,  the  fact  that  identical  simu¬ 
lated  pressure  fields  can  be  produced  for  a  range  of  resolution  conditions 
has  enabled  this  investigation  to  examine  the  influence  of  spatial  resolu¬ 
tion  (or  "microphone"  si2e)  on  the  statistics.  That  is,  an  identical  pres¬ 
sure  history  can  be  used  to  produce  the  pressure  force  history  for  different 
sized  receivers,  thereby  providing  an  unambiguous  measure  of  the  influence 
of  sensor  size  on  the  statistics  that  are  measured. 

Schewe  (Ref.  10)  has  measured  turbulent  wall  pressure  fluctuations 
using  an  array  of  transducers  of  varying  size.  His  measurements  suggest 
that  as  the  resolution  area  is  increased  the  statistical  properties  of  the 
local  pressure  force  approach  Gaussian.  Much  of  the  pressure  correlation 
data  has  been  summarized  by  Willmarth  (Ref.  11)  and  Von  Winkle  et  al.  (Ref. 
6)  and  the  reader  is  referred  to  those  reviews  for  additional  references. 
In  the  discussion  which  follows,  only  those  experimental  measurements  which 
were  used  in  developing  or  testing  the  simulation  will  be  cited.  In  addi¬ 
tion,  the  previous  work  of  the  authors  (Refs.  12,  13)  may  be  helpful  in 
developing  a  better  understanding  of  the  approach. 
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2.  SIMULATION  STRATEGY 


The  simulation  program  was  developed  as  a  multipurpose  program  to  model 
a  variety  of  flow  and  geometry  conditions.  The  program  is  listed  in 
Khorrami  (Ref.  13)  and  is  made  flexible  by  assuming  that  pressure  fluctua- 
;  tion  statistics  scale  according  to  the  parameters  listed  in  Table  1.  That 

table  also  contains  the  flow  variables  used  in  the  example  simulations.  The 
simulated  flow  conditions  were  similar  to  the  experiments  of  Bull  (Ref.  14), 

I  but  several  other  flow  conditions  have  been  examined  subsequently  (Ref.  2, 

15). 

The  simulation  record  is  constructed  from  randomly  distriouted,  dis- 
j  Crete  pressure  events.  Single  cycle  pressure  fluctuations  were  generated 

using  empirically  based  distribution  functions  for  spacing,  time  of  origin, 
peak  amplitude,  and  wavelength.  Depending  upon  the  frequency  (wavelength), 

I  the  individual  events  were  assigned  convection  speeds  and  convected  down¬ 

stream  across  the  user  specified  measurement  locations.  The  force  resulting 
from  the  event  was  computed  over  the  specified  measurement  area  by  integra- 
I  tion  at  the  specified  sampling  times.  The  instantaneous  pressure  force  was 

I 

the  SLTimation  of  all  the  events  that  contributed  to  that  measurement  area  at 
the  prescribed  time. 

I  A  symbolic  representation  of  the  pressure  events  that  are  created  dur¬ 

ing  a  single  spatial  pass  down  a  streamwise  track  is  shown  in  Fig.  1.  Each 
Ax^  is  generated  randomly  using  a  distribution  function  developed 
I  empirically  from  the  burst  cycle  data  of  Of fen  and  Kline  (Ref.  16).  The 

individual  event  location  is  determined  as  the  sun  of  the  Ax^'s  that  are 
generated  during  that  sweep.  When  the  location  of  origin  exceeds  the  total 
I  length  of  the  numerical  experiment,  a  new  sweep  is  started.  The  incremental 

time  of  origin  is  also  based  on  the  data  of  Offen  and  Kline  (Ref.  16),  but 


the  event  generation  times  were  related  to  a  reference  or  datum  time,  rather 
than  accumulated  during  the  sweep.  At  the  end  of  a  spatial  sweep,  an  aver¬ 
age  time  interval  is  calculated  for  that  sweep  and  the  reference  or  datum 
time  is  increased  by  the  average  time  step,  before  the  next  sweep  is 
started.  That  method  prevented  any  inadvertant  frequency  biasing  due  to 
unnecessary  time  step  control.  Obviously,  the  user  specified  resolution 
time  is  not  connected  to  the  simulation  time  and,  in  effect,  a  freely  vary¬ 
ing  clock  and  a  control  clock  are  used  simultaneously  to  manage  tne  super¬ 
position  calculations. 

Besides  the  known  variation  of  convection  speed  with  pressure  fluctu¬ 
ation  frequency,  experimental  evidence  suggests  that  the  frequency  or  wave¬ 
length  of  the  event  also  varies  as  the  fluctuation  is  convected  downstream 
(Refs,  6,  14,  and  17).  A  wavelength  adjustment  was  made  continuously  as  the 
event  convected  down  stream. 

A  characteristic  convection  speed  or  celerity  can  be  determined  by 
dividing  the  separation  distance  between  two  measurement  points  by  the  time 
interval  associated  with  the  peak  in  the  two-point  correlation  function.  It 
was  observed  during  the  development  phase  that  simulations  produced  convec¬ 
tion  speeds  which  exceeded  free  stream  velocity  when  adjacent  points  were 

★ 

close  together  <  2],  even  though  the  assigned  event  con¬ 
vection  speeds  were  always  less  than  0.8U^.  That  effect  was  attributable 
to  generation  of  pressure  events  which  spanned  more  than  one  measurement 
location  at  the  time  of  origin.  When  new  events  actually  bridged  two 
locations,  an  "infinite"  convection  speed  contribution  could  be  produced. 
That  problem  was  eliminated  by  ramping  the  amplituo-.  of  the  event  up  to  the 
prescribed  value  during  a  short  time  interval  [which  was  characteristic  of 
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the  time  associated  with  a  bursting  event  in  the  Offen  and  Kline  (Ref.  16) 
study).  Subsequently,  the  amplitude  of  the  pressure  disturbance  began  to 
decay.  A  schematic  representation  of  the  evolution  and  convection  of  a 
simulated  pressure  event  is  shown  in  Fig.  2. 

Because  of  the  non-linear  character  of  both  the  simulation  events  and 
the  statistical  behavior  of  the  turbulent  pressure  field,  the  simultaneous 
representation  of  pressure  statistics  (r.m.s.  pressure,  power  spectrum  and 
two  point  correlation)  involved  a  great  deal  of  tedious  trial  and  error  ad¬ 
justments  of  the  distribution  functions.  An  acceptable  combination  of  dis¬ 
tribution  functions  was  evolved  and  they  are  contained  in  the  computer  pro¬ 
gram  which  is  listed  in  Reference  (13).  A  flow  chart  describing  the 
simulation  program  is  shown  in  Figure  3. 

The  simulation  which  has  just  been  described  is  designed  to  generate  an 

arbitrarily  resolved  record,  P(x.  t  ),  in  the  streamwise  direction.  Since 

T ,  n 

other  records  over  the  same  spatial  interval  and  time  can  be  created  simply 
by  changing  the  initial  random  nunber,  a  set  of  pressure  records  --  say 
P^(x^.  t^)  --  can  be  produced.  Depending  on  the  length  of  the  time  record 
and  the  randomness  of  the  random  nunber  generator,  it  is  possible  to  con¬ 
struct  a  set  of  pressure  records  for  spatial  locations,  x^ ,  which  are 
statistically  independent.  That  is. 


N-k 

^(x,t)P  (x,t  +t)  =  I  P  (x.  t  )  P  (x.  t„  +  kAt)/[(N  -  k)P„  P.  >0  (1) 

^  n=l^  "  s  j,  n  s^^^ 

for  all  combinations  of  r,  s,  i,  j  and  k,  except  r  =  s.  Here,  N  is  the 
total  number  of  time  steps,  and  we  require  that  0  <  k  <  N/2.  As  N  is  in¬ 
creased,  the  correlations  should  approach  zero.  The  ability  to  produce  a 
set  of  statistically  independent  simulation  records  is  required  to  construct 
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a  complete  P(x,  z,  t)  simulation. 


A.  Simulation  of  P(x,  z,  t) 


If  the  pressure  field  is  to  be  resolved  at  uniformly  spaced  locations, 

z.  in  the  spanwise  direction,  there  is  some  distance,  L  ,  in  which  the 
J  J  ^ 

correlation  between  spanwise  pressure  records  is  important.  From  the  stand¬ 
point  of  the  simulation,  the  specified  length  of  the  time  record  controls 
the  deviation  of  the  correlations  in  Eq.  (1).  A  numerical  scatter  in  the 
correlations  is  present  for  finite  records  and  that  scatter  is  used  as  a 
basis  for  delineating  the  correlation,  width,  L^.  That  is,  if 


PJx,  t)  P(k,  t  +  t)  <  e 


for  time  records  of  length  NAt,  then  there  is  a  corresponding  interval, 
L^,  such  that 


P(x,  z,  t)  P(x,  z  +  Az  ,  t  +  t)  >  e. 


for  AZjjj  less  than  or  equal  to  L^. 

The  experimental  data  of  Bull  (Ref.  14)  and  Willmarth  and  Wooldridge 
(Ref.  17)  can  be  used  to  estimate  the  interval  over  which  a  spanwise 
correlation  is  greater  than  or  equal  to  e. 

For  a  spanwise  spacing  of  Az,  there  are  statistically  connected 
pressure  simulation  tracks. 


^n^’  ’*n  +  i^’^i*  ^n^*  ‘“•^n  +  ^n^’ 


where  is  the  integer  value  of 

This  control  requirement  becomes  quite  involved  for  a  rectangular  area  which 
is  wide  compared  to  Az,  since  the  correlations  must  be  controlled  for  all 
adjacent  simulation  tracks  in  both  lateral  directions.  The  correlations  can 
be  generated  by  constructing  P(x^,  Zj,  t^)  from  more  than  one  simulation 
track.  That  is,  we  can  let 

Nz 

P(«,.  Zj.  .  i(x,.  »„).  (5) 


where 


Nz 


I 

M=1 


1, 


(6) 


and 


ai  32  +  az  33  +  ...  +  3^^  _  j 


+  3, 


n. 


Ri 


ai  83  +  82  3^  +  ...  + 


Rz 

R 


N,  -  1, 


(7) 


The  values  of  R^ ,  R2 _ R,ij  .  1  3i"e  specified  from  the  experimental  span- 

★ 

wise  correlation  data  (Ref.  14,  17)  for  separation  distances  of  az/6  , 
2a z/ 6  ,  etc.  The  values  of  the  coefficients,  a^^,  must  be  determined  by 
trial  and  error.  It  follows  that 
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P(x.,  Zi,  t^)  =  ai  Pi(x..  tj  +  32  P2(x.,  t^)  +  ...  +  3^  P^  (Xi,  t^) 


PCX^,  Z2,  =  3i  P2(X.,  t^)  +  32  P2(X^.  t^)  + 


An  example  simulation  has  been  described  previously  (Ref.  12)  where  a 


spanwise  separation  distance  of  2.5  5  was  used  and  the  time  record  (4000 


data  points  per  station)  was  such  that  a  scatter  of  approximately  10  percent 


was  observed  in  the  uncontrolled  correlation  functions  (e  =  0.10).  Under 


*  * 

those  conditions  the  correlations  at  2.5  6  and  5  6  were  controlled,  while 


correlations  for  larger  distances  were  uncontrolled  and,  therefore,  scat¬ 


tered  about  ±0.10.  Variation  in  correlations  from  the  simulation  are  shown 


in  Fig.  4,  along  with  a  nominal  fit  of  the  experimental  data  for  t  =  0 


(Ref.  14,  17). 


3.  RESULTS  AND  DISCUSSION 

A  typical  segment  of  simulated  pressure  is  shown  in  Fig.  5  for  three 
locations  over  the  same  time  interval.  For  small  separation  distances 
(ax^  =  770  or  Ax/6*  =  1.7),  the  pressure  histories  are  similar,  but  at 
Ax^=  6100,  little  similarity  remains.  The  receiver  "diameter"  was  u^  d/v  = 
36  or  d/6  =  0.08^.  Schewe  (Ref.  10)  has  shown  similar  data  for  actual 
turbulent  wall  pressure  measurements. 

It  is  also  interesting  to  observe  the  influence  of  the  receiver  area  on 
the  recorded  signal  driven  by  an  identical  pressure  field.  The  pressure 
histories  generated  on  a  receiver  area  (at  x=o)  which  is  half  as  large  as 
the  one  in  Figure  5  (i.e.,  d^  =  18)  and  for  one  which  is  more  than  an  order 
of  magnitude  larger  (d^  =  270)  are  shown  in  Figure  6.  The  reference  case 
(d^  =  36)  is  also  shown.  While  the  records  for  the  two  small  receivers  are 
very  similar,  it  is  obvious  that  the  large  receiver  area  results  in  an 
integrated  force  per  unit  area  record  which  is  very  different. 

The  characteristic  dimension  of  the  pressure  receivers  can  be  changed 
simply  by  changing  a  line  of  computer  code.  By  subjecting  different  sensing 
areas  to  the  same  pressure  history,  the  influence  of  receiver  size  on  root 
mean  square  pressure  can  be  studied.  For  u^  d/v  >  400,  spanwise  decorre¬ 
lation  must  be  considered,  but  with  that  effect  neglected,  r.m.s.  pressure 
is  shown  as  a  function  of  receiver  size  in  Fig.  7.  Those  data  appear  to 
agree  with  the  measurements  of  Bull  and  Thomas  (Ref.  18)  in  terms  of  the 
increase  in  r.m.s.  pressure  with  decreasing  transducer  size.  However,  the 
diameter  effect  starts  at  a  much  smaller  value  of  u^d/v  than  any  of  the 
experimental  measurements. 

^For  a  characteristic  sensing  area  dimension,  d/6  <  1,  the  spanwise 
correlation  measurements  are  sufficiently  close  to  unity  to  ignore  spanwise 
effects.  Then  the  streamwise  length  is  equivalent  to  a  diameter. 
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Wall  pressure  energy  spectra  are  shown  for  three  different  receiver 
areas  in  Fig.  8  along  with  a  cross  hatched  composite  of  experimental  mea- 
suranents.  The  peak  frequency  becomes  delineated  more  clearly  as  the  sens¬ 
ing  area  becomes  smaller,  and  there  is  a  slight  increase  in  high  frequency 
energy.  However,  overall  shape  of  the  power  spectrum  is  not  sensitive  to 
receiver  area. 

Variation  of  streamwise  two  point  correlation  functions  with  receiver 
size  and  separation  distance  are  shown  in  Fig.  9.  It  is  interesting  to  note 
that  the  correlations  spread  and  diminish  as  the  receiver  size  is  decreased. 
This  effect  is  due  to  the  inclusion  of  small  scale  pressure  fluctuations 
which  decay  more  rapidly  and  thus  decorrelate  adjacent  pressure  signals.  A 
slight  irregularity  in  the  correlations  for  the  smallest  receiver  area  is 
due  to  the  prescribed  time  step  resolution  and  should  be  ignored. 

Because  of  the  way  in  which  the  simulation  was  produced  it  was  easy  to 
employ  a  pseudo-filter  and  examine  the  streamwise  correlations.  If  a  random 
event  had  a  frequency  of  interest,  it  was  allowed  to  contribute  to  the  wall 
pressure  data.  If  it  did  not  have  a  frequency  of  interest,  it  was  ignored. 
Thus,  two  point  correlations  could  be  developed  which  consisted  of  only 
those  events  containing  prescribed  frequencies.  High  and  low  frequency 
correlations  are  shown  for  a  range  of  separation  distances  in  Figs.  10  and 
11  along  with  the  experimental  measurements  of  Willmarth  and  Wooldridge 
(Ref.  17).  Since  single  cycle  events  carry  side  band  frequencies,  in  addi¬ 
tion  to  their  prescribed  frequencies,  this  filtering  technique  is  not  per¬ 
fect. 

Schewe  (Ref.  10)  examined  the  influence  of  receiver  size  on  the  skew¬ 
ness  and  flatness  of  the  wall  pressure  fluctuation  signal.  A  comparison 
between  his  measurements  and  the  simulation  data  are  shown  in  Figs.  12  and 


13.  While  the  skewness  data  agree  that  there  is  a  higher  probability  for 
negative  pressure  fluctuations  as  the  receiver  area  decreases,  there  is  a 
significant  discrepancy  between  the  flatness  or  kurtosis  of  the  simulation 
and  Schewe's  measurements. 


4.  CONCLUSIONS 


A  Monte  Carlo  based  simulation  has  been  developed  which  can  produce 
simulated  turbulent  wall  pressure  fields.  The  statistics  generated  from  the 
simulations  are  in  reasonably  good  agreement  with  experimental  measurements 
and  since  the  simulated  field  can  be  replicated  for  different  receiver  di¬ 
mensions,  it  was  possible  to  compare  statistical  properties  for  systems 
subjected  to  identical  pressure  fields. 

—  /  Based  on  these  simulations,  it  was  possible  to  conclude  that  the  root 
mean  square  pressure  levels  increase  in  a  quasi  linear  manner  as  the  re¬ 
ceiver  size  ("microphone")'  decreases.  The  trend  is  in  substantial  agreement 
with  the  experiments  of  Bull  and  Thomas,,  (Ref .  but  the  threshold  of  the 
diameter  effect  and  the  magnitude  of  the  r.m.s.  increase  may  be  controlled 
by  flow  phenomena  that  are  either  ignored  or  improperly  simulated.  The 
power  spectra  are  nearly  insensitive  to  receiver  size  in  the  energy  contain¬ 
ing  frequency  interval.  Two-point  correlations  first  show  higher  correla¬ 
tions  with  decreasing  receiver  size,  then  show  poorer  correlations  as  the 
receiver  size  becomes  small  enough  to  sense  fine  scale  phenomena. 

The  authors  believe  this  simulation^ program  can  be  a  valuable  tool^in 
studying  the  response  of  complex  or  non-linear  structures  to  quasi-random 
wall  pressure  fields.  The  ability  to  adjust  resolution  and  simulated  flow 
conditions  arbitrarily  makefjttrrs' a  flexible  tool  in  the  analys+s' and  design'^ 
ef' fluid-structural  systems. 
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Variation  of  the  amplitude  and  wavelength  of  a  typical 
simulated  pressure  event  as  it  is  convected  downstream 


3.  Flow  chart  representing  the  program  for  simulating  P(x,t). 
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risen  of  lateral  correlations  from  simulated  P(x,z,t)  and  the  experimental  data 
references  14  and  17. 


downstream  locations  (Ax 


Figure  8.  Influence  of  measurement  area  on  power  spectra  for 
identical  pressure  histories. 
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Comparison  of  filtered  two-po1nt  correlations  with  experimental  data 
High  frequency  case. 


Willmarth  and  Wooldridge  (17) 


Comparison  of  filtered  two- point  correl 
experimental  data.  Low  frequency  case. 


